Comparison of Flow Reduction Efficacy of Nominal and Oversized Flow Diverters Using a Novel Measurement-assisted in Silico Method.

PURPOSE: The high efficacy of flow diverters (FD) in the case of wide-neck aneurysms is well demonstrated, yet new challenges have arisen because of reported posttreatment failures and the growing number of new generation of devices. Our aim is to present a measurement-supported in silico workflow that automates the virtual deployment and subsequent hemodynamic analysis of FDs. In this work, the objective is to analyze the effects of FD deployment variability of two manufacturers on posttreatment flow reduction. METHODS: The virtual deployment procedure is based on detailed mechanical calibration of the flow diverters, while the flow representation is based on hydrodynamic resistance (HR) measurements. Computational fluid dynamic simulations resulted in 5 untreated and 80 virtually treated scenarios, including 2 FD designs in nominal and oversized deployment states. The simulated aneurysmal velocity reduction (AMVR) is correlated with the HR values and deployment scenarios. RESULTS: The linear HR coefficient and AMVR revealed a power-law relationship considering all 80 deployments. In nominal deployment scenarios, a significantly larger average AMVR was obtained (60.3%) for the 64-wire FDs than for 48-wire FDs (51.9%). In oversized deployments, the average AMVR was almost the same for 64-wire and 48-wire device types, 27.5% and 25.7%, respectively. CONCLUSION: The applicability of our numerical workflow was demonstrated, also in large-scale hemodynamic investigations. The study revealed a robust power-law relationship between a HR coefficient and AMVR. Furthermore, the 64 wire configurations in nominal sizing produced a significantly higher posttreatment flow reduction, replicating the results of other in vitro studies.

Fractals and Chaos in the Hemodynamics of Intracranial Aneurysms.

Computing the emerging flow in blood vessel sections by means of computational fluid dynamics is an often applied practice in hemodynamics research. One particular area for such investigations is related to the cerebral aneurysms, since their formation, pathogenesis, and the risk of a potential rupture may be flow-related. We present a study on the behavior of small advected particles in cerebral vessel sections in the presence of aneurysmal malformations. These malformations cause strong flow disturbances driving the system toward chaotic behavior. Within these flows, the particle trajectories can form a fractal structure, the properties of which are measurable by quantitative techniques. The measurable quantities are well established chaotic properties, such as the Lyapunov exponent, escape rate, and information dimension. Based on these findings, we propose that chaotic flow within blood vessels in the vicinity of the aneurysm might be relevant for the pathogenesis and development of this malformation.

Image-based flow simulation of platelet aggregates under different shear rates.

Hemodynamics is crucial for the activation and aggregation of platelets in response to flow-induced shear. In this paper, a novel image-based computational model simulating blood flow through and around platelet aggregates is presented. The microstructure of aggregates was captured by two different modalities of microscopy images of in vitro whole blood perfusion experiments in microfluidic chambers coated with collagen. One set of images captured the geometry of the aggregate outline, while the other employed platelet labelling to infer the internal density. The platelet aggregates were modelled as a porous medium, the permeability of which was calculated with the Kozeny-Carman equation. The computational model was subsequently applied to study hemodynamics inside and around the platelet aggregates. The blood flow velocity, shear stress and kinetic force exerted on the aggregates were investigated and compared under 800 s-1, 1600 s-1 and 4000 s-1 wall shear rates. The advection-diffusion balance of agonist transport inside the platelet aggregates was also evaluated by local Péclet number. The findings show that the transport of agonists is not only affected by the shear rate but also significantly influenced by the microstructure of the aggregates. Moreover, large kinetic forces were found at the transition zone from shell to core of the aggregates, which could contribute to identifying the boundary between the shell and the core. The shear rate and the rate of elongation flow were investigated as well. The results imply that the emerging shapes of aggregates are highly correlated to the shear rate and the rate of elongation. The framework provides a way to incorporate the internal microstructure of the aggregates into the computational model and yields a better understanding of the hemodynamics and physiology of platelet aggregates, hence laying the foundation for predicting aggregation and deformation under different flow conditions.

Modelling collateral flow and thrombus permeability during acute ischaemic stroke.

The presence of collaterals and high thrombus permeability are associated with good functional outcomes after an acute ischaemic stroke. We aim to understand the combined effect of the collaterals and thrombus permeability on cerebral blood flow during an acute ischaemic stroke. A cerebral blood flow model including the leptomeningeal collateral circulation is used to simulate cerebral blood flow during an acute ischaemic stroke. The collateral circulation is varied to capture the collateral scores: absent, poor, moderate and good. Measurements of the transit time, void fraction and thrombus length in acute ischaemic stroke patients are used to estimate thrombus permeability. Estimated thrombus permeability ranges between 10-7 and 10-4 mm2. Measured flow rates through the thrombus are small and the effect of a permeable thrombus on brain perfusion during stroke is small compared with the effect of collaterals. Our simulations suggest that the collaterals are a dominant factor in the resulting infarct volume after a stroke.

Loss of α4A- and ß1-tubulins leads to severe platelet spherocytosis and strongly impairs hemostasis in mice.

Native circulating blood platelets present with a discoid flat morphology maintained by a submembranous peripheral ring of microtubules, named marginal band. The functional importance of this particular shape is still debated, but it was initially hypothesized to facilitate platelet interaction with the injured vessel wall and to contribute to hemostasis. The importance of the platelet discoid morphology has since been questioned on the absence of clear bleeding tendency in mice lacking the platelet-specific ß1-tubulin isotype, which exhibits platelets with a thinner marginal band and an ovoid shape. Here, we generated a mouse model inactivated for ß1-tubulin and α4A-tubulin, an α-tubulin isotype strongly enriched in platelets. These mice present with fully spherical platelets completely devoid of a marginal band. In contrast to the single knockouts, the double deletion resulted in a severe bleeding defect in a tail-clipping assay, which was not corrected by increasing the platelet count to normal values by the thrombopoietin-analog romiplostim. In vivo, thrombus formation was almost abolished in a ferric chloride-injury model, with only a thin layer of loosely packed platelets, and mice were protected against death in a model of thromboembolism. In vitro, platelets adhered less efficiently and formed smaller-sized and loosely assembled aggregates when perfused over von Willebrand factor and collagen matrices. In conclusion, this study shows that blood platelets require 2 unique α- and ß-tubulin isotypes to acquire their characteristic discoid morphology. Lack of these 2 isotypes has a deleterious effect on flow-dependent aggregate formation and stability, leading to a severe bleeding disorder.

Traumatic vessel injuries initiating hemostasis generate high shear conditions.

Blood flow is a major regulator of hemostasis and arterial thrombosis. The current view is that low and intermediate flows occur in intact healthy vessels, whereas high shear levels (>2000 s-1) are reached in stenosed arteries, notably during thrombosis. To date, the shear rates occurring at the edge of a lesion in an otherwise healthy vessel are nevertheless unknown. The aim of this work was to measure the shear rates prevailing in wounds in a context relevant to hemostasis. Three models of vessel puncture and transection were developed and characterized for a study that was implemented in mice and humans. Doppler probe measurements supplemented by a computational model revealed that shear rates at the edge of a wound reached high values, with medians of 22 000 s-1, 25 000 s-1, and 7000 s-1 after puncture of the murine carotid artery, aorta, or saphenous vein, respectively. Similar shear levels were observed after transection of the mouse spermatic artery. These results were confirmed in a human venous puncture model, where shear rates in a catheter implanted in the cubital vein reached 2000 to 27 000 s-1. In all models, the high shear conditions were accompanied by elevated levels of elongational flow exceeding 1000 s-1. In the puncture model, the shear rates decreased steeply with increasing injury size. This phenomenon could be explained by the low hydrodynamic resistance of the injuries as compared with that of the downstream vessel network. These findings show that high shear rates (>3000 s-1) are relevant to hemostasis and not exclusive to arterial thrombosis.

The effect of stiffened diabetic red blood cells on wall shear stress in a reconstructed 3D microaneurysm.

Blood flow within the vasculature of the retina has been found to influence the progression of diabetic retinopathy. In this research cell resolved blood flow simulations are used to study the pulsatile flow of whole blood through a segmented retinal microaneurysm. Images were collected using adaptive optics optical coherence tomography of the retina of a patient with diabetic retinopathy, and a sidewall (sacciform) microaneurysm was segmented from the volumetric data. The original microaneurysm neck width was varied to produce two additional aneurysm geometries in order to probe the influence of neck width on the transport of red blood cells and platelets into the aneurysm. Red blood cell membrane stiffness was also increased to resolve the impact of rigid red blood cells, as a result of diabetes, in blood flow. Wall shear stress and wall shear stress gradients were calculated throughout the aneurysm domains, and the quantification of the influence of the red blood cells is presented. Average wall shear stress and wall shear stress gradients increased due to the increase of red blood cell membrane stiffness. Stiffened red blood cells were also found to induce higher local wall shear stress and wall shear stress gradients as they passed through the leading and draining parental vessels. Stiffened red blood cells were found to penetrate the aneurysm sac more than healthy red blood cells, as well as decreasing the margination of platelets to the vessel walls of the parental vessel, which caused a decrease in platelet penetration into the aneurysm sac.

The Effects of Micro-vessel Curvature Induced Elongational Flows on Platelet Adhesion.

The emerging profile of blood flow and the cross-sectional distribution of blood cells have far reaching biological consequences in various diseases and vital internal processes, such as platelet adhesion. The effects of several essential blood flow parameters, such as red blood cell free layer width, wall shear rate, and hematocrit on platelet adhesion were previously explored to great lengths in straight geometries. In the current work, the effects of channel curvature on cellular blood flow are investigated by simulating the accurate cellular movement and interaction of red blood cells and platelets in a half-arc channel for multiple wall shear rate and hematocrit values. The results show significant differences in the emerging shear rate values and distributions between the inner and outer arc of the channel curve, while the cell distributions remain predominantly uninfluenced. The simulation predictions are also compared to experimental platelet adhesion in a similar curved geometry. The inner side of the arc shows elevated platelet adhesion intensity at high wall shear rate, which correlates with increased shear rate and shear rate gradient sites in the simulation. Furthermore, since the platelet availability for binding seems uninfluenced by the curvature, these effects might influence the binding mechanics rather than the probability. The presence of elongational flows is detected in the simulations and the link to increased platelet adhesion is discussed in the experimental results.

Modelling the leptomeningeal collateral circulation during acute ischaemic stroke.

A novel model of the leptomeningeal collateral circulation is created by combining data from multiple sources with statistical scaling laws. The extent of the collateral circulation is varied by defining a collateral vessel probability. Blood flow and pressure are simulated using a one-dimensional steady state blood flow model. The leptomeningeal collateral vessels provide significant flow during a stroke. The pressure drop over an occlusion predicted by the model ranges between 60 and 85 mmHg depending on the extent of the collateral circulation. The linear transport of contrast material was simulated in the circulatory network. The time delay of peak contrast over an occlusion is 3.3 s in the model, and 2.1 s (IQR 0.8-4.0 s) when measured in dynamic CTA data of acute ischaemic stroke patients. Modelling the leptomeningeal collateral circulation could lead to better estimates of infarct volume and patient outcome.

Biorheology of occlusive thrombi formation under high shear: in vitro growth and shrinkage.

Occlusive thrombi formed under high flow shear rates develop very rapidly in arteries and may lead to myocardial infarction or stroke. Rapid platelet accumulation (RPA) and occlusion of platelet-rich thrombi and clot shrinkage have been studied after flow arrest. However, the influence of margination and shear rate on occlusive clot formation is not fully understood yet. In this study, the influence of flow on the growth and shrinkage of a clot is investigated. Whole blood (WB) and platelet-rich plasma (PRP) were perfused at high shear rates (> 3,000 s-1) through two microfluidic systems with a stenotic section under constant pressure. The stenotic section of the two devices are different in stenotic length (1,000 vs 150 µm) and contraction angle of the stenosis (15° vs 80°). In all experiments, the flow chamber occluded in the stenotic section. Besides a significantly increased lag time and decreased RPA rate for PRP compared to WB (p < 0.01), the device with a shorter stenotic section and steeper contraction angle showed a shear-dependent occlusion and lag time for both PRP and WB. This shear-dependent behavior of the platelet aggregate formation might be caused by the stenotic geometry.

A novel virtual flow diverter implantation method with realistic deployment mechanics and validated force response.

Virtual flow diverter deployment techniques underwent significant development during the last couple of years. Each existing technique displays advantageous features, as well as significant limitations. One common drawback is the lack of quantitative validation of the mechanics of the device. In the following work, we present a new spring-mass-based method with validated mechanical responses that combines many of the useful capabilities of previous techniques. The structure of the virtual braids naturally incorporates the axial length changes as a function of the local expansion diameter. The force response of the model was calibrated using the measured response of real FDs. The mechanics of the model allows to replicate the expansion process during deployment, including additional effects such as the push-pull technique that is required for the deployment of braided FDs to achieve full opening and proper wall apposition. Furthermore, it is a computationally highly efficient solution that requires little pre-processing and has a run-time of a few seconds on a general laptop and thus allows for exploratory analyses. The model was applied in a patient-specific geometry, where corresponding accurate control measurements in a 3D-printed model were also available. The analysis shows the effects of FD oversizing and push-pull application on the radial expansion, surface density, and on the wall contact pressure.

The influence of red blood cell deformability on hematocrit profiles and platelet margination.

The influence of red blood cell (RBC) deformability in whole blood on platelet margination is investigated using confocal microscopy measurements of flowing human blood and cell resolved blood flow simulations. Fluorescent platelet concentrations at the wall of a glass chamber are measured using confocal microscopy with flowing human blood containing varying healthy-to-stiff RBC fractions. A decrease is observed in the fluorescent platelet signal at the wall due to the increase of stiffened RBCs in flow, suggesting a decrease of platelet margination due to an increased fraction of stiffened RBCs present in the flow. In order to resolve the influence of stiffened RBCs on platelet concentration at the channel wall, cell-pair and bulk flow simulations are performed. For homogeneous collisions between RBC pairs, a decrease in final displacement after a collision with increasing membrane stiffness is observed. In heterogeneous collisions between healthy and stiff RBC pairs, it is found that the stiffened RBC is displaced most. The influence of RBC deformability on collisions between RBCs and platelets was found to be negligible due to their size and mass difference. For a straight vessel geometry with varying healthy-to-stiff RBC ratios, a decrease was observed in the red blood cell-free layer and platelet margination due to an increase in stiffened RBCs present in flow.

Hydrodynamic Resistance of Intracranial Flow-Diverter Stents: Measurement Description and Data Evaluation.

PURPOSE: Intracranial aneurysms are malformations forming bulges on the walls of brain arteries. A flow diverter device is a fine braided wire structure used for the endovascular treatment of brain aneurysms. This work presents a rig and a protocol for the measurement of the hydrodynamic resistance of flow diverter stents. Hydrodynamic resistance is interpreted here as the pressure loss versus volumetric flow rate function through the mesh structure. The difficulty of the measurement is the very low flow rate range and the extreme sensitivity to contamination and disturbances. METHODS: Rigorous attention was paid to reproducibility, hence a strict protocol was designed to ensure controlled circumstances and accuracy. Somewhat unusually, the history of the development of the rig, including the pitfalls was included in the paper. In addition to the hydrodynamic resistance measurements, the geometrical properties-metallic surface area, pore density, deployed and unconstrained length and diameter-of the stent deployment were measured. RESULTS: Based on our evaluation method a confidence band can be determined for a given deployment scenario. Collectively analysing the hydrodynamic resistance and the geometric indices, a deeper understanding of an implantation can be obtained. Our results suggest that to correctly interpret the hydrodynamic resistance of a scenario, the deployment length has to be considered. To demonstrate the applicability of the measurement, as a pilot study the results of four intracranial flow diverter stents of two types and sizes have been reported in this work. The results of these measurements even on this small sample size provide valuable information on differences between stent types and deployment scenarios.

Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH)-phase II: rupture risk assessment.

PURPOSE: Assessing the rupture probability of intracranial aneurysms (IAs) remains challenging. Therefore, hemodynamic simulations are increasingly applied toward supporting physicians during treatment planning. However, due to several assumptions, the clinical acceptance of these methods remains limited. METHODS: To provide an overview of state-of-the-art blood flow simulation capabilities, the Multiple Aneurysms AnaTomy CHallenge 2018 (MATCH) was conducted. Seventeen research groups from all over the world performed segmentations and hemodynamic simulations to identify the ruptured aneurysm in a patient harboring five IAs. Although simulation setups revealed good similarity, clear differences exist with respect to the analysis of aneurysm shape and blood flow results. Most groups (12/71%) included morphological and hemodynamic parameters in their analysis, with aspect ratio and wall shear stress as the most popular candidates, respectively. RESULTS: The majority of groups (7/41%) selected the largest aneurysm as being the ruptured one. Four (24%) of the participating groups were able to correctly select the ruptured aneurysm, while three groups (18%) ranked the ruptured aneurysm as the second most probable. Successful selections were based on the integration of clinically relevant information such as the aneurysm site, as well as advanced rupture probability models considering multiple parameters. Additionally, flow characteristics such as the quantification of inflow jets and the identification of multiple vortices led to correct predictions. CONCLUSIONS: MATCH compares state-of-the-art image-based blood flow simulation approaches to assess the rupture risk of IAs. Furthermore, this challenge highlights the importance of multivariate analyses by combining clinically relevant metadata with advanced morphological and hemodynamic quantification.

A new hypothesis on the role of vessel topology in cerebral aneurysm initiation.

Aneurysm pathogenesis is thought to be strongly linked with hemodynamical effects. According to our current knowledge, the formation process is initiated by locally disturbed flow conditions. The aim of the current work is to provide a numerical investigation on the role of the flow field at the stage of the initiation, before the aneurysm formation. Digitally reconstructed pre-aneurysmal geometries are used to examine correlations of the flow patterns to the location and direction of the aneurysms formed later. We argue that a very specific rotational flow pattern is present in all the investigated cases marking the location of the later aneurysm and that these flow patterns provide the mechanical load on the wall that can lead to a destructive remodelling in the vessel wall. Furthermore, these patterns induce elevated vessel surface related variables (e.g. wall shear stress (WSS), wall shear stress gradient (WSSG) and oscillatory shear index (OSI)), in agreement with the previous findings. We emphasise that the analysis of the flow patterns provides a deeper insight and a more robust numerical methodology compared to the sole examination of the aforementioned surface quantities.

Cellular Level In-silico Modeling of Blood Rheology with An Improved Material Model for Red Blood Cells.

Many of the intriguing properties of blood originate from its cellular nature. Therefore, accurate modeling of blood flow related phenomena requires a description of the dynamics at the level of individual cells. This, however, presents several computational challenges that can only be addressed by high performance computing. We present Hemocell, a parallel computing framework which implements validated mechanical models for red blood cells and is capable of reproducing the emergent transport characteristics of such a complex cellular system. It is computationally capable of handling large domain sizes, thus it is able to bridge the cell-based micro-scale and macroscopic domains. We introduce a new material model for resolving the mechanical responses of red blood cell membranes under various flow conditions and compare it with a well established model. Our new constitutive model has similar accuracy under relaxed flow conditions, however, it performs better for shear rates over 1,500 s-1. We also introduce a new method to generate randomized initial conditions for dense mixtures of different cell types free of initial positioning artifacts.

Towards the virtual artery: a multiscale model for vascular physiology at the physics-chemistry-biology interface.

This discussion paper introduces the concept of the Virtual Artery as a multiscale model for arterial physiology and pathologies at the physics-chemistry-biology (PCB) interface. The cellular level is identified as the mesoscopic level, and we argue that by coupling cell-based models with other relevant models on the macro- and microscale, a versatile model of arterial health and disease can be composed. We review the necessary ingredients, both models of arteries at many different scales, as well as generic methods to compose multiscale models. Next, we discuss how this can be combined into the virtual artery. Finally, we argue that the concept of models at the PCB interface could or perhaps should become a powerful paradigm, not only as in our case for studying physiology, but also for many other systems that have such PCB interfaces.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Emerging fractal patterns in a real 3D cerebral aneurysm.

The behaviour of biological fluid flows is often investigated in medical practice to draw conclusions on the physiological or pathological conditions of the considered organs. One area where such investigations are proven to be useful is the flow-related formation and growth of different pathologic malformations of the cerebro-vascular system. In this work, a detailed study is presented on the effect of a cerebral aneurysm on blood transport inside a human brain artery segment. This malformation causes strong flow instabilities that drives the flow system towards chaotic behaviour. The emerging fractal structure and some of its measurable properties have been explored using a method that makes the measurement of these properties feasible even in complicated large three dimensional data sets. We find that, from the investigated chaos parameters, the information dimension turns out to be the most reliable parameter to characterise chaotic advection in the vicinity of the aneurysm sac. We propose that properties of chaotic mixing close to aneurysms might be relevant for the condition of this pathologic malformation.